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The inter-band optical absorption in graphene characterized by its Hne-structure constant has a 
universal value of 2.3% independent of the material parameters. However, for several graphene- 
based photonic applications, enhanced optical absorption in graphene is highly desired. In this 
work, we quantify the tunability of optical absorption in graphene via the Fermi level in graphene, 
angle of incidence of the incident polarized light, and the dielectric constant of the surrounding 
dielectric media in which graphene is embedded. The influence of impurities adsorbed on the 
surface of graphene on the Lorentzian broadening of the spectral function of the density of states is 
analytically evaluated within the equilibrium Green’s function formalism. Finally, we compute the 
differential absorption of right and left circularly-polarized light in graphene that is uniaxially and 
optically strained. The preferential absorption or circular dichroism is investigated for armchair and 
zigzag strain. 


I. Introduction 

The two-dimensional material graphen^i^^, which is 
a layer of carbon atoms arranged in a hone yco mb lat¬ 
tice, exhibits strong light-matter interactiorP^ over a 
very wide wavelength ranging from the far infrared to 
the ultraviolet. The tunability of the density of states 
and the Fermi level in graphene along with its excellent 
transport properties reflected in a high carrier mobilitj^^ 
provide a pat h for photonic app lications such as quan¬ 
tum optics^, photo- voltaicd^lio)^ photo-detectoriS^, and 
biological sensing.E^ While the optical characteristics of 
mono-layer graphene, attributed to enhanced light mat¬ 
ter absorption with a high quantum efficiency, make it 
a desirable material for optical resonators and thermal- 
imaging cameras, the o ptical absorption is poor to be an 
efficient photo-detector.li^Illl Further, the adaptability of 
the optical absorption in graphene for various frequencies 
is limited by a flat abso:^ion spectrum in the visible to 
the near-infrared region.lISI 

In this work, we propose methods to improve the over¬ 
all absorption in mono-layer graphene in a dielectric en¬ 
vironment via a direct tuning of its optical conductivity. 
We begin with a description of the optical characteris¬ 
tics of a graphene sheet suspended in vacuum and derive 
the quantum of light absorption for a beam at oblique 
incidence (Section II A). The calculations are repeated 
for a graphene sheet sandwiched between two dielectrics; 
the dielectric function of graphene in each case is estab¬ 
lished using a standard calculation. While most 

calculations tacitly assume an idealized set-up in which 
graphene is pristine thus preserving its electronic struc¬ 
ture, especially the linearly dispersing bands around the 
K and K edges of the Brilloui n zon e, impurity atoms are 
usually adsorbed on the surfac4^^^*^to alter the electronic 
spectrum and attendant optical response. The correction 
in this situation is manifest through a spectral broaden¬ 
ing of the density of states and a changed absorption pro¬ 
file. The broadening for an impurity-adsorbed graphene 
sheet is evaluated by employing a Hubbard-type Hamil¬ 


tonian within the equilibrium Green’s function formalism 
(Section H B). Finally, in Section H C, a uniaxially- and 
optically-strained graphene sheet is analyzed for circular 
dichroism. Results are collected in Section HI, and the 
paper concludes by briefly touching upon the implica¬ 
tions of increased absorption for graphene-based optical 
devices. 


II. Optical absorption in suspended graphene 


Dispersion relationships for graphene bands with lin¬ 
early dispersing eigen stateJi^ are modeled using a two- 
dimensional Dirac Hamiltonian given as 


PL = hvj{axky - ayk,^) -I- Acr^, (1) 

where Vf denotes the Fermi velocity of carriers, k is 
the wave-vector measured relative to the Dirac points, 
and Ui {i = x, y, z) are the usual Pauli matrices. The 
corresponding wave functions, in momentum space, for 
the momentum around the Dirac point have the form: 


where rj 
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/ Urj{k) exp{—i9)\ 
rju-yik) ) ’ 


±, and Uy is given by 


(2a) 


Uy{k) 



A 

y^A^ -I- (hvfk)^ 


(2b) 


Here A is the band gap induced in graphene, 9k = 


arctan 



and ?7 = ± corresponds to conduction 


and valence bands. The band gap introduced is primar¬ 
ily through an interaction with the substrate on which 
graphene is eptaxially grown, for instance, th e hon ey- 
comb lattice rigidly held on the hexagonal ex¬ 

hibits a band gap between 7.0 meV to 20.0 meV.I^ 
The band gap of pristine graphene, which hosts mass¬ 
less Dirac fermions, is zero. 
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Graphene-dielectric interface 

FIG. 1: Schematic of a TM (p-polarized) wave incident at 
angle di on the graphene-dielectric interface. The incident 
ray is partly reflected and partly transmitted. The angles of 
reflection and transmittance are indicated on the sketch. The 
electric field is normal to the propagation vector, while the 
magnetic field is along the t/-axis. It is assumed that there is 
no beam incident on the interface from medium 2. 


A. Transfer matrix for optical absorption 

To determine the optical response of carriers in 
graphene, linearly-polarized light along x-axis is assumed 
to shine perpendicularly on the graphene surface as in¬ 
dicated in Fig. The electric field of the light beam 
can be determined using the vector potential ~^{t) = 
exp^—iujt), where uj is the frequency of the incident 
light. The corresponding electric field is given as 




( 3 ) 


The speed of light is denoted by c in Eq. The 
Hamiltonian of Eq. can therefore be modified us¬ 
ing the Peierls substitutioiP^ and takes the form H — 
hvf~^. (^k — We extract the interaction part of 

the Hamiltonian from the modified Hamiltonian and is 
expressed as 


T~Lint = — 


hvfC 




( 4 ) 

2oj 


Substituting for from Eq. 3 T-Lmt = 

The factor of 0.5 in Eq. comes from by retaining only 
the real part in the expansion of the magnetic vector 
potential 

Using Fermi golden rule, the transition probability for 
light-induced transition from valence to conduction band 
is 

fcc.fct, 


energies at the bottom of the conduction band and the 
top of the valence band, respectively. 

Recognizing the delta function in Eq. as the density of 
2|e| 

states given as . „ „ for the linear energy-dispersion re- 
lationship in graphene and integrating over energy space 

g 2|^|2 

gives the absorbed flux (see Appendix A) as ———. The 

incident flujP^ for an electric field is —\E\^. The opti- 

4_7r 

cal inter-band absorption in graphene is, therefore, uni- 
2 

7r6^ 

versally given as = 2.3%, a material-independent 
nc 

number. 

The optical absorption changes when graphene is sand¬ 
wiched between two dielectric layers characterized by di¬ 
electric constants of Ci and 62 - A transverse-magneticP^ 
(TM or p - polarized wave) with a magnetic field along 
the j/-axis (Eq. is assumed to impinge on the graphene- 
dielectric interface. The magnetic field corresponding to 
the TM polarization and propagating along the x-axis is 
given as 






( 6 ) 


where kz = kin cos 0 and kx = km sin 0. A particular 
physical quantity within the «*^layer in the structure is 
identified by the subscript i, km is incident wave vector, 
and the angle of incidence is 0. The appropriate bound¬ 
ary conditions at the interface z = 0 are 


and 


n X (l^j —l^i^ =0, 

^ j ^ 2 ^ — surf 7 


(7a) 


(7b) 


where surf is the surface current due to the 2D charge 
carriers in the graphene sheet. The subscripts i, j denote 
a particular material layer on the left and the right side 
of a given interface, respectively. The electric field can 

^ 

be obtained using Maxwell’s equation V xE = —— to 

at 


yield 


= - —^xt. 


iO£o£i 


For the form of the magnetic 


field chosen in Eq. the corresponding electric field is 
given as 


Ez,. = 


kx 


.z,^ - -(8a) 


weoG 


Ex,i — 


ojeod 




—ikz iZ 


D ikz ■iZ\ 

-B,e 


(8b) 


Combining the above equations with the boundary con¬ 
dition (Eq. 7a) for electric field on either side of the 
graphene-dielectric interface gives 


Here, ipf and ipi correspond to the final and the initial 
scattering states, respectively; Ec and Ey refer to the 


fcz.l 

ei 


(Ai 


Bi) 


kz ,2 


(A2 


B2) = 0 . 


(9a) 


£2 
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The second relation using boundary conditions applied 
on the magnetic field (Eq. 7b) is 


Jx — (^1 + Bx) — (A 2 + B 2 ) , 
kz,2 (^2 ~ B 2 ) 

= <7 -. 

£0^2 


(9b) 

(9c) 


Note that in the above equation, tr is the complex dy¬ 
namical conductivity of graphene as discussed in Section 
II B. All quantities on left(right) of the interface are sub¬ 
scripted as “1(2)”. In matrix notation, the amplitude 
of the magnetic fields on either side of the interface are 
related ad^ 



( 10 a) 


Since all quantities in the above equation are assumed to 
be positive, the left half side of the above inequality is 
trivially true; evaluating the right hand side inequality. 


£2 


^2 


one arrives at the relation — < —. For a case which 

_£i fci , 

is in violation of this condition, unphysical solutions are 
obtained. We also demonstrate a parallel condition by 
placing a similar constraint on the reflectance amplitude. 
Proceeding as above, we write 


0 < r = 


Mtm (2,1) 

Mt™(l,l) 


< 1 


(13c) 


Ignoring the contribution of Mt m ,2 i n the high-frequency 
limit, the right inequality {Eq. 13c I is trivially satisfied 
while the left takes the form 


The transfer matrix Mtm = Mi^tm + 1^2,trm where 


Ml 



(l + 

£2^z,1 

_ £lfcz,2 
V f-2kz,l 


1 - 
1 + 


eikz,2 \ 

£2^z,1 

eikz,2 

e2kz,i J 


and 


M 2 



/ 0'fcz,2 

£o£2W 

crkz,2 

\£o£2W 


O-fcz.2 

£o£2W 

CTfcz,2 

£o£2W 


( 10 b) 


(lOc) 


The transfer matrix M^^ derived in Eq. |10a| allows us 
to compute reflectance and transmittance amplitudes; in 
particular, reflectance and transmittance amplitudes are 


( 2 , 1 ) 

1 

Mi™ (1,1)' 


( 11 ) 

( 12 ) 


The reflection and transmission coefficients are R = 
and T = Y which add to unity (i? -I- T =1) for zero 
absorption losses. 

At this point it is imperative to discuss appropriate 
limits that are placed on the expressions for reflectance 
and transmittance amplitudes derived in Eqs. [TT]and[T^ 
The amplitudes must satisfy 0 < r < 1 and 0 < t < 1 
from which follows 


0 < t 


1 

Mi™ (1,1) 


< 1 . 


(13a) 


To evaluate the constraint, we note that the matrix 
Mtm ,2 which describes the optical conductivity of the 
graphene sheet makes a small contribution to the overall 
matrix Mtm in the high-frequency limit. We are, there¬ 
fore, left with the inequality 

0 < t < 1. (13b) 

£2ftlz + £lfc2z 


£2fclz — £lfc2z 
£2fclz + f-lk2z ' 


(13d) 


It is straightforward to see that the condition r > 0 leads 
us to — > —, which means the transmittance ampli- 

ti ki 

tude is greater than unity. This apparent contradiction 

is resolved by noting that in this limiting case — = 7 —, 

, , £1 ki 

the reflected component is zero and the incident beam 
is fully transmitted {t = 1). We have thus arrived at 
an analogous condition for obtaining the corresponding 
Brewster’ angle for a graphene sheet embedded in inho¬ 
mogeneous dielectric media. Of course, the limits on the 
validity of the ratio € 2/^1 must also account for the finite 
contribution from the transfer matrix M 2 ,tm in the lower 
frequency regime. 

The transfer matrix Mte for a teP wave or s - polar¬ 
ized wave can be similarly derived using the appropri¬ 
ate Maxwell’s boundary conditions. The matrix Mte = 
Ml te + M 2 te is 



1 

' 2 

^J-lkz,2 




Mi^te - 

IJ-2kz,l 
, ^J‘lkz,2 

A< 2 fcz,l 
j 1 k-lkz,2 

5 

(14a) 


\ ^J-2kz,l 

M 2 fcz,l / 



and 







1 

" 2 


/Xo/riwcr \ 



M2,te 

^l.Z 

fioPi^cr 

Jolllwa 


(14b) 



V ^ 1 .^ 

ki,. J 



The magnetic permeabilities 

(/ri,/r 2 ) are 

taken 

to be 


unity for non-magnetic dielectric media. 


B. Optical absorption in impure graphene 


Electromagnetic absorption in graphene is determined 
by its complex dynamical conductivity. The dynamical 
conductivity of graphene is obtained from the dielectric 
function e {q, uj) through a random phase approximation 
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(RPA). In the long-wavelength limit {q 
dielectric function is given 


: (<? 0, w) = 1 - 


0), the RPA 


27re^ 

- 1 

1 


2Ef — Huj 

q 2TThui 

hoj 

h^ln 

2Ef + Huj 


-i-Q{hw-2Ef) 

(15) 

where Ef is the Fermi level in the graphene sheet, and 
0(.) is the Heaviside function. The dynamical conduc¬ 
tivity, a{q,uj), is related to the dielectric constant as 


ct(<7,w) = ^ [1 - e(9,w)]. 


(16) 


Inserting Eq. in Eq. gives 


a{q^0,u}) = 




le" , 
+ -.—y In 
AttH 


2Ef — huj 


2Ef -|- hui 


(17) 


The spectral width P corresp onds to a Lorentzian broad¬ 
ening of the density of state^^aiso] jg determined by 
the impurities adsorbed on the graphene surface. In this 
work, using a retarded Green’s function approach, we 
evaluate the spectral width, P. A model Hubbard Hamil- 
toniaiP^ that captures an atom adsorbed on graphene 
can be written as 




^ ^ ^ ^ ^hyh^. 


t 




r-E 


(18) 


hyb^s^gr- 


Expanding Eq.j^by evaluating the commutator (see Ap¬ 
pendix B), one obtains 


(c,s,0 = 6(t) + ^G^ {c,s,t)+^^VhybG^ {k,t), 


dt 


h 


( 21 ) 

where G^ {k,t) is the retarded Green’s function for 
graphene-adsorbed atom. G^ (fc, t) can be written in 
standard form as 

G^ {k, t) = -id (t) ({ogr (t), (0)}). (22) 

The corresponding equation of motion is 

(fc, t) = ^-fG^ (fc, t) + ^VhybG^ (c, s, t ). (23) 


The Fourier transform (Eq. [M[ ) of the two retarded 
Green’s functions yield a pair equations (Eq. 25 26) 
which can be solved for G^ (c, w) 

G^ic,s,t) = —j e-“*G«(c,s,a;), (24) 

{oj-^-fG^ (c, a;)) = 1 + ^Ky.G^ {k, oj ), (25) 


and 


(uj - ^G^ {k,ui)j = ^VhybG^ {c,s,u}). (26) 

Solving for G^ (c, Lo) and preserving causality by making 
the substitution uj uj + gives 


This Hamiltonian resembles the non-interacting Ander¬ 
son impurity model for resonant impurities.^^The energy 
of the adsorbed atom is given by Sad, while Sgr denotes 
the energy of the graphene electron states. The last two 
terms of the Hamiltonian describe the hybridization be¬ 
tween the adsorbed atom and graphene. The summation 
over momentum vectors also denote the dual spin states. 
The creation (annihilation) operators for graphene and 
the adsorbed atom are denoted a^(a) and cj(cs), respec¬ 
tively. 

To evaluate the spectral density function from which 
the overall broadening is determined, we consider the 
graphene-adsorbed impurity to be a non-interacting sys¬ 
tem. Using the equation of motion approach within 
the equilibrium Green’s function formalism, we start by 
writing the retarded Green’s functiorP^ for the adsorbed 
atom 

G^ (c, s, t) = -id {t) ({cs (t), (t)}). (19) 

Taking the derivative of the retarded Green’s func¬ 
tion in Eq. and following Heisenberg’s picture 
{H (t) = exp {iTit) % exp {—iHt)) gives 

^JG^ (c, s,t) = 6 {t) + ^0 (t) m, c. it )], 4 it)}). (20) 


G^ (c, s,a;) = 


huj + i0+ - Sad - E T— ■n-^ 
h 

huj -ead-J2 -1 *7r E \Vhyb\'^S {huJ - £gr) 

tlOJ — £ 


gr 


(27) 


where P stands for the principal valu^^ in the usual 

Plemelj relation: - = P ( — ) ± iirS {x). The delta 

X ±i£ \x J 

function is again the usual density of states for graphene; 
putting all of them together and evaluating the spectral 
density A (c, s, oj) = —2ImG^ (c, s, oj) gives 


A(c,s,w) = -- 


2;ir 


Hu^-Sad-EP^^ 

k nio — £g 


-, (28) 


■ r2 


where d{Ef) 


nlPv'j 


is the density of states for 


graphene at the Fermi level and the broadening pa¬ 
rameter F = TrS\Vhyb\‘^d{Ef); S denotes the area of 
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the graphene sheet under consideration and is roughly 
given as 5” « ^Iriimp^ where riimp is the concentra¬ 
tion of impurity atoms adsorbed on the graphene sheet. 
The hybridization potential can be estimated as Vhyb ~ 

hvf/ ln(i?o/i?i)|j.^Here, i?o is the radius of the 

impurity atom, and i?i is the average distance between 
the adsorbed impurity atoms and is approximately given 

^S 1/ yjTliYap- 

Impurities therefore broaden the Dirac-delta peak, the 
broadening given by the strength of the hybridization 
potential, Vhyb- We have tacitly assumed that the sys¬ 
tem is non-interacting, and there is no Coulomb repulsion 
term of the form that appears in a standard 

Hubbard Hamiltonian. 


where a± = (a^ ± icry) /2 and Eop is the optical strain. 
A gap, A = 2Eop, in the spectrum appears now which 
implies the non-equivalence of the graphene sub-lattices. 
Interestingly, the out-of-plane optical strain component, 
while it breaks the reflection symmetry, does not move 
the Dirac points from the Brillouin zone edges K and K . 

The optical matrix elements and the tunable circular 
dichroism of an optically strained graphene sheet can be 
calculated by first working out the velocity components. 
The X and y components of velocity are given by 


and 


dT-Lstr _ f 0 ^Vf\ 

~ 0 J ’ 


(33a) 


C. Circular dichroism in strained graphene 

Circular dichroism quantifies the differential absorp¬ 
tion of right and left circularly-polarized light. The ex¬ 
pression for absorption coefficient due to light-induced 
inter-band transitions is a meas ure of the strength of 
the optical matrix element Inter-band optical ab¬ 
sorption from a populated valence band eigen state to 
an empty conduction band state requires the determina¬ 
tion of the optical inter-band transition matrix elements 
|P+P and |Pci;P for right and left circularly-polarized 
light,respectively. for instance, for right circularly 

polarized light is nominally defined a^^ 

= (vI-+|Vfe,H|d/_). (29) 


dT-Latr _ / 0 —ihvf\ 
dky ~ yi^vf 0 ) ■ 


(33b) 


The optical matrix element corresponding to the velocity 
components for a finite band gap graphene can therefore 
be written by inserting expressions for wave functions 
and the velocity components from Eq. i and Eq. [3^ re¬ 
spectively, in Eq. [^for a right circularly-polarized light. 

^ 1 Kexp(«) „.) (4 -f') (“-") ■ 

(34) 

The matrix element corresponding to momentum opera¬ 
tor py is 


The degree of circular polarizatiorP^ p (fc) is therefore 


p{k) = 


IntP - \Pcv? 

\pXV + \pEvV' 


(30) 


The inter-band matrix elements for right and left 
circularly-polarized light, |P+| and |Pc«l) respectively are 
defined as 


P± = + ipy 

CV -^CV CD 


(31) 


We utilize Eq. 29 to compute the inter-band optical 
matrix elements in strained graphene and establish a re¬ 
lation between dischroism and the tunable Fermi level. A 
graphene sheet usually experiences acoustic and optical 
strain. Under acoustic strain, the two carbon atoms of 
the unit cell are displaced together, while for the optical 
case, there is a shift such that the centre-of-mass remains 
invariant. An out-of-plane component of optical strain 
renders the two sub-lattices of graphene inequivalent, a 
situation commonly realized when graphene is grown on 
a substrate, for instance, boron nitride. The Hamilto¬ 
nian using the method of invariants that describes this 
situation is written a^^ 


T-Latr = flVf {{kx - iky) Cr+ -I- {kx; + iky) tT_} -I- PopCT^ 


= kvf 


E. 


kx 'tk'i 

-^op 


op 

+ iky — E^ 


(32) 


(35) 

Combining both the components, the square of the right 
polarized optical matrix element is 


\Pcv? = {kV}f 


1 + ' 


A 


n 2 


(hvfk) + 


(36) 


The corresponding expression for following an 

analogous procedure is 


PcX = {kvfY 


1 + 


-A 


(hvfk) -I-A2 


(37) 


The above derivation for dichroism in graphene was 
carried out by considering an optical strain that furnishes 
a finite band gap. First principle calculation^^ also show 
that an application of a large uniaxial strain in graphene 
does not destroy the semi-metallic nature of graphene 
but significantly impacts the Fermi velocity components. 
The anisotropy of Fermi velocity which manifests as tilted 
Dirac cones has been profitably employed to tune the op¬ 
tical properties of two-dimensional graphene.^^ We work 
out an expression for circular dichroism in graphene by 
considering uniaxial strain along the zigzag and armchair 
directions^ in presence of a finite band gap. Writing 
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out the Hamiltonian (Eq. 321 again but with different 
velocity components, % = VxPxO'x + VyPyO'y, such that 
K = Vyjvx- The ratio is usually computed by resorting to 
a single-orbital tight-binding calculation which expresses 
the velocity components as^ 


Vx = haxV^v"^ - 1 ; 

Vy=t2ay. (38) 


Here, ti {i = 1,2,3) denotes the kinetic energy hopping 
integrals between the nearest neighbours, and ry = ti/t 2 - 
The strained lattice vector is resolved in x and y com- 

, -x/So 

ponents to give and Ux = a/2 and Uy = , with 


a = -s/S X I. 42 A. Using first principles calculationPi^, 
it has been demonstrated that a uniaxial strain along 
the zig-zag direction (denoted as “Z” strain) in the hon¬ 
eycomb lattice leads to ti = < t 2 { p < 1), while for 

uniaxial strain in the armchair chain direction (denoted 
as “A” strain), ti = > t 2 [p > 1). 

The degree of velocity anisotropy, k, using Eq. is 
given as 




(dry^ — 1) 


(39) 


III. Results 

As a first check of the validity of the numerical model 
to compute absorption, an undoped graphene sheet with 
an electrically tuned Fermi level and suspended in vac¬ 
uum (the dielectric constants flanking the graphene sheet 
are unity) is considered. The absorption coefficient, at 
normal angle of incidence, for a range of energies is plot¬ 
ted in Fig. Two features in Fig. stand out; when 
incident radiation energy is upwards of twice the Fermi 
level, the absorption is constant at approximately 2.3 % 
and the absorption coefficient has a “hump” at exactly 
hu) = 2Ef. The numerical demonstration of absorption 
coefficient as approximately 2.3 % is significant since it 
is in conformity with the theoretically derived number 
(see Appendix A) obtained directly from the wave func¬ 
tions of graphene’s linear Hamiltonian coupled to the in¬ 
cident electromagnetic field. The “hump” is attributed 
to an onset of inter-band conductivity at huj = 2Ef, 
beyond which the absorption coefficient stays constant. 
The intra-band Drude contribution is dominant at much 
lower energies as can be seen by observing the optical 
conductivity versus energy plot of Fig. The “tuning” 
of absorption coefficient via the optical conductivity of 
graphene will concern us for rest of this section. 


The degree of circular polarization (Eq. 301 is modified to 
reflect this by an alteration to the algebraic expressions 
for both right and left polarized optical matrix elements. 
Carrying out the derivation as before gives 




1-b 


±A 


^J {hVxkxf' + {hVykyf + A2_ 


(40) 

where Vr = \J'^x + terms of the velocity 

anisotropy, the above equation can be simplified as 


1 + 


= (f^Vx) {1 + K^)x 

-I 2 

±A 


{hvxkxf (1 -I- tan^ 9k) + A^ 


(41) 


where tan Ok = ky/kx- 

The average degree of circular polarization over a con¬ 
stant energy surface (free of trigonal warping effect^^Slieij 
can be defined as 



FIG. 2: The absorption coefficient of mono-layer graphene 
suspended in vacuum is plotted for various energies. Here, 
ujf IS the frequency corresponding to the Fermi level in the 
graphene sheet. As shown, a constant 2.3% absorption is 
obtained independent of the material parameters for uj > ujf. 
While the Fermi level, Ef, plays a role in determining absorp¬ 
tion in the low-frequency regime where a; << a;/, in the high- 
frequency regime characterized by inter-band scatterings, ab¬ 
sorption coefficient becomes independent ot Ef. 


^ J d^kp{k)S{uj -2eik)) 

’ Jd^k6{oj-2e{k)) 


(42) A. Absorption coefficient versus incident angle and 

material parameters 


where e (fc) = EyJWv/k^PpA^ is the energy spectrum of 
graphene. The frequency of the light beam must satisfy 
Lo = 2e (k) which describes the energy involved in an 
inter-band transition from valence to conduction band. 


The amount of light absorbed by a suspended graphene 
sheet must be a tunable quantity for a wide range of 
applications such as photo-detectors and sensors. This 
tuning, for a given set of material constants, can be ac¬ 
complished by selecting “control” parameters such as fre- 
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FIG. 3: Complex dynamical conductivity of mono-layer 
graphene. The Fermi level is set to 0.2 eV. Fig. plots 
the Drude (intra-band) component of conductivity, while the 
inter-band component is shown in Fig. [2>. The damping fac¬ 
tor is assumed to be 2.6 meV. The real (imaginary) part of 
each contribution is a solid (dashed) line in each sub-figure. 
The a;-axis is normalized to frequency corresponding to the 
Fermi level. The y-axis is normalized to cto = 7re^/2/i. 



FIG. 4: Absorption coefficient is shown as a function of the 
incident angle. Different traces correspond to the marked 
ratio ^ = ei/e 2 of the dielectric constants. Other simulation 
parameters are noted in the figure legend. 


quency of the incoming light-beam, the Fermi level and 
the angle of incidence to alter the overall reflectance, 
transmission, and absorption. We show in Fig. the 
relation between incident angle, 9, and the absorption 
coefficient for various values of ^ = ei/e 2 at fixed val¬ 
ues of Ef^ F, and w. The two dielectric constants sur¬ 
rounding the graphene sheet are denoted by ei and €2 as 
sketched in Fig. It is easily seen that the absorption 
coefficient reaches up to 35.0 % when the ratio of dielec¬ 
tric constants is 2.0 at close to normal incidence. Besides 
the enhanced absorption, another noteworthy feature in 
Fig. ID is the degradation in absorption coefficient as the 
incident angle increases from zero to 7r/2. This is ex¬ 
plained by examining the expressions for reflectance and 
transmittance given in Eqs. |13b[ |13d[ for a fixed ratio of 
the dielectric constants, at higher angles of incidence, the 
transmitted wave vector (in medium “2”) ^ 2,2 changes 
such that a greater portion of the incident light is re¬ 
flected. The transmitted wave in medium “2” is simply 

expressed as ^ 2,2 = (^ — 'j £2 — where = ko sin 6 . 

ko and w describe the incident wave-vector and frequency, 
respectively. A lesser absorption with decreasing ratio of 

the dielectric constants — is also explained by utilizing 
£2 

Eqs. |13b[ |13d[ it is evident that a lower ratio as marked 
on the plot will augment the reflection coefficient. 

We next turn our attention to absorption characteris¬ 
tics as a function of the Fermi level and the operating 
frequency. The functional dependence is shown in Fig. 
[^for absorption at normal incidence, £1 = 4, £2 = 2, and 
the broadening parameter F set to 2.6 meV. At an oper¬ 
ating frequency of around ui « O.Olw/, the absorption co¬ 
efficient exhibits a peak for Ef — 0.4 eV and Ef = 0.6 eV 
as shown in Fig.[^ Two crucial observations can be made 
that will help us define an “optimal” parameter space to 
design graphene photonic devices. We first focus on the 


low bjjbjj ratio which gives an absorption peak; at this 
operational frequency, the Drude inter-band conductiv¬ 
ity dominates such that absorption increases up to ap¬ 
proximately 70.0 %, saturates, and then starts to fall. 
For w > O.lujf, the Drude conductivity which is still the 
dominant mode until w = 2a;/, rapidly drops to make a 
negligible contribution independent of the Fermi level(see 
Fig.i; with all other parameters held constant, we are 
therefore able to explain the merging of the absorption 
profiles for all three Fermi levels. The inter-band scat¬ 
tering is zero since we increase a; until it reaches a;/. 

The inset in Fig. [^further shows the role of the broad¬ 
ening parameter in determining the absorption of light. 
As the broadening parameter gains strength, in regions 
where Drude conductivity is significant (a; < O.lo;/), the 
conductivity is lowered for identical operational charac¬ 
teristics such as Fermi level, frequency etc. The lowered 
conductivity explains the drop in absorption coefficient. 
The broadening parameter in the plot is a Lorentzian and 
is assumed to be Fermi level independent. 

We have thus identified a small window for the oper¬ 
ating frequency for Ef > 0.4 eV, which enables a max¬ 
imized absorption up to approximately 70.0 % for the 
selected dielectric constants surrounding the graphene 
sheet for TM polarized incident light. A low value of 
broadening parameter is desirable to keep the absorption 
coefficient high. While in the above estimates for absorp¬ 
tion, the broadening was an empirically chosen number, 
we present below a more accurate calculation that re¬ 
lates it to the impurity concentration and Fermi level of 
graphene. 

An impure graphene sheet, which is a more realistic 
scenario, is considered to evaluate the spectral broaden¬ 
ing of states. The impurities are non-interacting as de¬ 
scribed by the Hamiltonian in Eq. [T^ Figure conveys 
the role of impurities in influencing the optical absorp¬ 
tion, where the Lorentzian broadening is characterized by 




















FIG. 5: Absorption coefficient is shown as a function of the 
operating frequency normalized to the frequency correspond¬ 
ing to the Fermi level in graphene. Different traces correspond 
to the marked values of Fermi level in graphene. The inset 
plot shows the impact of Lorentzian broadening on the ab¬ 
sorption coefficient in graphene. 


the hybridization potential, of the impurity atom 

with the TT-bonds of graphene. The trend in Fig. is in 
conformity with the behaviour shown in the inset plot 
of Fig. As expected, with an increase in the impurity 
concentration, the broadening is enhanced which lowers 
the absorption. The three trend lines shown in Fig. [^fur¬ 
ther demonstrate that for a low w/w/ ratio, the Drude 
conductivity is significant as indicated by the downward 
trend in absorption(a;/a;/ = 0.05). As w —>■ w/, Drude 
conductivity is considerably decreased in magnitude as 
depicted by the “relatively” straight absorption curves. 
These absorption curves are also insensitive to the Fermi 
level in graphene. However, in the presence of other im¬ 
purity atoms that are interacting such as charged impuri¬ 
ties and polar phonons, absorption coefficient will indeed 
exhibit a dependence on the Fermi level. This aspect has 
not been considered. The hybridization potential used 
in determining the broadening is shown in the inset plot 
of Fig. [^as a function of riimp- For a typical impurity 
concentration of riimp ~ 10^° Vhyb ~ (0.5 — 1) 

KeVA^P 

We have therefore identified a few key guidelines to 
design a graphene photonic device. Primarily, in the 
low-frequency regime where Drude conductivity is ac¬ 
tive, absorption saturates as a function of the ratio of 
dielectric constants and frequency of incident light. The 
corresponding Fermi level can be adjusted such that ab¬ 
sorption is maximized. Finally, the broadening parame¬ 
ter must not be too high, a large spectral broadening at 
low frequencies degrades the overall absorption. 

B. Circular Dichroism 

Circular dichroism or the degree of circular polariza¬ 
tion relates to the differential absorption of right and 
left circularly-polarized light. In the case of graphene 



FIG. 6: Absorption coefficient versus impurity concentration 
in impure graphene for different values of operating frequency 
normalized to the frequency corresponding to Fermi level. 
The inset plot shows the hybridization potential modeled as 
a function of impurity concentration. 

with its linear bands, the dichroism is analytically com¬ 
puted through Eq. As is clear from Eq. pristine 
graphene (A = 0) does not exhibit circular polarization 
p\ however, in the presence of an inversion symmetry¬ 
breaking band gap, which renders the two sub-lattices 
inequivalent, p is a non-zero number. We seek to evalu¬ 
ate p in the following section under a symmetry-breaking 
condition. The coupling of right (left) circularly polar¬ 
ized light (see Fig. which is exact at the Dirac point 

formed at K (^K ^ edge of the Brillouin zon^^ exhibits 

a polarization-dependent light absorption near a particu¬ 
lar edge, say K. This differential absorption quantified as 
the degree of circular polarization is plotted against the 
Fermi level. Th e Fermi level in this case is assumed to 
coincide with ^(hvfk)^ -\- A^. In principle, though the 
degree of circular polarization is a function of the mo¬ 
mentum vector, the variation of p is shown for a single k 
point in Fig. 




FIG. 7: Right and left circularly polarized light couples selec¬ 
tively to band edges K and K since they are time-reversed 
pairs. At other points in fc-space, not far from the band 
edge, the selective coupling is lost and a varying degree of 
polarization-dependent absorption occurs that gives rise to 
the phenomenon of circular dichroism. The two curved ar¬ 
rows denote left and right circular polarization. 

There are several instances where graphene is not only 
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FIG. 8: Degree of circular polarization p in gapped graphene 
versus the Fermi level for two values of the band-gap: A = 5 
meV and 10 meV. p is plotted using Eq. over a constant 
energy surface dehned by the respective Fermi level, which is 
dehned by the energy spectrum of graphene. The fc-vector is 
taken as 0.1 1/A. The dichroism decreases with an increasing 
Fermi level. 


band gap open (A ^ Q) K and K edges but also asym¬ 
metrically strained (uniaxial strain) to produce tilted 
Dirac cones. To understand circular dichroism in this sit¬ 
uation, we first plot in Fig. [fusing data from Ref.lHl the 
anisotropy of x— and y— directed components of Fermi 
velocity of the Dirac fermions in a uniaxially “A” strained 
graphene. The marked difference in the velocity compo¬ 
nents is evident, with an increase in “A” strain, the y- 
component of the Fermi velocity decreases, while the x- 
component is boosted such that the velocity-anisotropy 
factor K = Vyjvx reduces to 0.35 from unity. The in¬ 
crease in uniaxial strain until 24% does not contribute to 
an additional band gap other than that created through 
substrate-induced optical strain. Choi et al. in Ref. 01 ] 
explain this reduction in k by considering the strength 
of the hopping integral parameters. The qualitative rela¬ 
tionship between hopping integrals, connecting the near¬ 
est neighbors in a single-orbital tight-binding approxima¬ 
tion under “A” strain in graphene, ti = > t 2 is shown 

in the inset of Fig. Using the velocity-anisotropy graph 
of Fig. the degree of circular polarization in uniaxially- 
strained graphene is computed (Fig. 101 as a function of 
“A” strain percentage. We immediately notice from the 

k 

plot that for 9k = 7r/4, where 9k = arctan the degree 

kx 

of circular polarization does not considerably vary while 
there is a noticeable change for 9 = tt/O, 7 r/3 with increas¬ 
ing strain. The trend for d = 7r/4 is explainable if we note 
that at this angle, the two components of the A:-vectors 
are symmetric, the uniaxial strain does not contribute 
to the band gap A and the velocity anisotropy-factor 
only causes a slight change in degree of circular polar¬ 
ization. The slight deviation as strain increases and con¬ 
sequently the velocity-anisotropy factor k supports this 
reasoning. The other angles introduce additional asym¬ 


metry that augments the velocity-anisotropy to produce 
a larger change in the degree of circular polarization. 

The same behaviour is observed again for strain of the 
zig-zag type which is plotted in Fig. 11 At 0 = 7r/4, the 
degree of circular polarization is almost constant while 
it shows an upward trend for 9 — tt/6 . The important 
point to note is the inter-play of /c-space and velocity 
anisotropy; if the two degrees of asymmetry cancel out, 
as it happens for 0 = 7r/3 in the zig-zag case, there will 
not be a sufficient variation in circular poalrization. This 
assumes importance in view of our proposal to enhance 
light absorption of a particular polarization, by select¬ 
ing a suitable A:-space asymmetry combined with velocity 
anisotropy we can tune the differential absorption of right 
or left circularly-polarized light. This idea runs parallel 
to the emergin g field of valleytronics in transition metal 

dichalcogenidesPSHni 



FIG. 9: Data on velocities Vy and Vx versus uniaxial “A” 
strain as extracted from calculations in Ref. SB The inset 
shows the honeycomb lattice of graphene under strain, where 
the kinetic energy hopping integrals ti are shown. In this case, 
ti = t 2 > ts. The inset schematically depicts distortion of the 
graphene honeycomb lattice under uniaxial armchair strain. 


IV. Conclusions 

We have demonstrated in this work the tunability of 
optical absorption in graphene via a combination of fab¬ 
rication design techniques and external dynamic control. 
The insertion of dielectric films of specified permitivity 
that surround the active graphene layer govern the trans¬ 
mission and reflection coefficients at a given Fermi level. 
The reflection and transmission coefficients which can 
then be made to span the complete range, R,T £ [0,1], 
therefore allow, through a selection of the dielectric con¬ 
stants, a desired absorption window. This window can 
be adjusted to conform to any sought level of absorption 
for a particular application. We refer to the choice of di¬ 
electric layers as a possible fabrication design technique. 
The effect of surface adsorbed impurity atoms is also a 
viable mechanism to tune the optica l conductivity; by a 
careful choice of impurity atoms^^^, the spectral broad- 



















10 



FIG. 10: Degree of circular polarization in graphene ver¬ 
sus the uniaxial armchair strain for two different values of 
bandgap: A = 5 meV and 10 meV. The degree of circular 
polarization depends on the anisotropy of k space and Fermi 
velocity components. For a larger band gap, the differential 
absorption of light is signihcant as borne out by the two well- 
separated group of lines for A = 5,10 meV. 


graphene. We have not considered in this work the influ¬ 
ence of strain that changes the atomic orbital overlap to 
give rise to an anisotropic conductivity tensor. Strain was 
considered in an elementary treatment to measure the de¬ 
gree of circular polarization assuming that the Dirac dis¬ 
persion is preserved; in reality, strain alters the density 
of states, the Dirac Hamiltonian and underlying optical 
and electronic properties. A more complete investigation 
of optimally engineered strained graphene structures will 
be considered in a follow-up work. Further, controllable 
doping in graphene, which provides a pathway for easy 
switching^ between n-type and p-type can be considered 
as a fabrication design parameter to enhance optical ab¬ 
sorption. Doping has not been explicitly considered in 
this work. 
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FIG. 11: Degree of circular polarization in graphene versus 
the uniaxial zigzag strain for two different values of bandgap: 
A = 5 meV and A = 10 meV. For 6 = tt/S, the degree of 
circular polarization does not show any variation. 


ening of the density of states can be adjusted to match 
the proper optical absorption. Circular dichroism is con¬ 
sidered as another possible option to selectively absorb 
polarization-dependent light by focusing attention to op¬ 
tical processes in one of the two edges (^K,K'^. Real 
time dynamic control, different from fabrication designs 
can also be exercised by a gate contact that alters the 
Fermi level in the graphene sheet. A changing Fermi 
level influences graphene’s optical conductivity and can 
switch the graphene-based optical device between vary¬ 
ing degrees of absorption via modification of the reflec¬ 
tion and transmission coefficients. 

There are several parameters that can be explored 
to further refine the efficiency of optical absorption in 


Appendix A Optical absorption in mono-layer 

graphene 

For the sake of completeness, we derive the absorbance 
for a graphene sheet and show that it is independent 
of the incident photon frequency. The incident light is 
assumed to be linearly polarized along x-axis and shines 
normally on the surface. The electric field, using the 
vector potential ^(t) = ~^exp{—iujt) is 




\dl 

c dt 


(Al) 


The modified Hamiltonian after including the Peierls sub¬ 
stitution takes the form % = Vfl^ ■ The 

interaction part of the Hamiltonian is therefore Hint = 

2 c 

Substituting for ^ from Eq. 


Al 


n. 


= The 

2tLl 


factor of 0.5 comes from by retaining only the {—iujt) 
term. 

Using the Fermi golden rule, the transition probability 

for a carrier to be excited from the valence band to con- 
Stt 

duction band is X! f\'^int\'^i)\‘^5 {Ec - Ey - hw) 


where the delta function is transformed in to the density 
of states. 

In the above expressions, iki and represent the initial 
and final wave functions while E^. and Ey are the energies 
corresponding to the bottom of the conduction band and 
top of the valence band, respectively. Calculating the 
matrix element M (k) = f\'Hint\'^i) gives 


M(/c) = i(exp(10) . (A2) 
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Inserting the expression for T-Lint in Eq. X2 and noting 
that the Ux dots with the x polarized electric field, one 
obtains 


M{k) 


Vfe\Ex 

Aiuj 

Vfe\Ex 

2ui 


(exp(i0) 
sin 6*, 



^exp(-i6»)^ 

(A3) 


where the wave functions of Eq. are utilized. The 
probability density by inserting the square of the matrix 
element in Fermi Golden rule can be written as 


Ph-fe 


2nvy\E^\ 

h 


sin^ 9 


(fr. 


(A4) 


The absorbed energy flux density is therefore (j) = 
— hwPh^e- Summing over the density of states which is 


2e 


, where energy e = 


huj 

~2 


and using the mean value 


of sind over a full cycle as 1/2, the final energy flux takes 
the following expression 


This Ti. gives four commutators. The first commutator is 




^ (a^ {ay,cj- |a^,Cs|ay) . 

(B2) 


Expanding terms within the bracket, the commutator is 
zero. The relation \AB , C\ = A{B,C} — {A,C}B is 
used here. In Eq. B2 the creation (annihilation) op¬ 
erators for graphene and the adsorbed atom are simply 
denoted by aj (og) and (cg) respectively. The subscript 
s describes the spin projection. More detailed subscripts, 
on-site energies of graphene and adsorbed impurity atom, 
and hybridization potentials that appear in the original 
Hamiltonian (Eq. B11 are omitted for brevity. 

The second commutator is 


E = E (^!' - {^I' ’ Cg'). 


(B3) 


1 27r 1 hw 

(ptotal = 4^2 2 2TTh?v‘] 

_ 

“ Ah ■ 


(A5) 

(A6) 


Note that the integral f (Pr is equal to A, the exposed 
graphene area in Eq. |A4| The incident flux for an electric 

C 7 r 6 ^ 

field I if I is —|ifp- The absorption is therefore = 
47r he 

2.3%. This number includes the valley degeneracy of two 
in graphene. 


Appendix B 


Evaluation of the commutator in 
Eq. 


20 


The final expression for the commutator in Eq. 
is derived here. The commutator is given as 
{[if, Cg (t)], c| (0)}. We proceed by simplifying the in¬ 
ner commutator [i?, Cg (t)]. The Hamiltonian if is given 
by 


T~L — ^ ^QT^gr^gr T ^ ^ ( ^hyb^, 




A, 


(BI) 


hyb^s^gr- 


Standard commutation relations in Eq. |B4| simplify the 
commutator^ c^c^',Cg to —Cg . 


{c^SCg} = 0; {cl,Cg} = I. 


(B4) 


The third commutator from the Hamiltonian is 






Evaluating each sub-anticommutator on RHS and using 
commutation relations in Eq. |B4| yields zero. 

The last commutator is 






As before, the evaluating the sub-anticommutators give 


c\, a^> , Cg 
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